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I. INTRODUCCION 



El metodo canonico de cuantizacion permite establecer una correspondencia entre las 
cantidades clasicas observables y los operadores cuanticos. El hecho de que al emplearlo 
para cuantizar campos aparezcan energias de vacio divergentes es una consecuencia de la 
ambigiiedad en la eleccion del orden usado para escribir los terminos del liamiltoniano, que 
son productos de operadores que no conmutan entre si. Los hamiltonianos escritos con 
diferentes ordenes difieren en terminos proporcionales al operador identidad Restar 
estos terminos para pasar de un ordenamiento a otro es equivalente a trasladar el nivel cero 
de energia. Ya que solamente las diferencias de energias entre estados son observables, la 
eleccion del cero de energia no afecta la dinamica de los campos. Asi, en la teoria de campos 
libres, la ambigiiedad se resuelve con la prescripcion de orden normal que conduce a un valor 
medio de vacio nulo para el operador liamiltoniano. 

En 1948, H. B. G. Casimir mostro que la presencia de dos placas conductoras neutras 
conduce a un valor medio de vacio no nulo para el operador liamiltoniano de un campo 
electromagnetico 

Por extension, se conoce como efecto Casimir a la variacion que sufren los valores medios 
de vacio de observables de un campo al imponersele condiciones de contorno externas o al 
interactuar con un campo externo generado por determinada configuracion de fuentes. La 
diferencia entre las energias de vacio de estas dos configuraciones se define como Energia de 
Casimir. Asi, se tiene 

Ecas = Eq — Eq. (1) 

donde Eq es la energia de vacio del campo libre y Eo es la energia de vacio del campo 
distorsionada por condiciones de contorno o la interaccion con un campo externo. 

Los metodos de evaluacion de energias de vacio pueden reunirse en dos categorias 
0: el metodo de suma directa de modos, aplicable toda vez que el espectro clasico sea 
explicitamente determinable, y el metodo local, que requiere el conocimiento de la funcion 
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de Green del problema. En este trabajo, usaremos el metodo de suma de modos. Segiin este, 
la energia de vacio se obtiene de sumar los autovalores de energia de los infinites modos del 
problema particular. Para el caso de un campo de Dirac como el que trataremos, la energia 
de vacfo resulta |Q 

EQ = -\Y,gk\Eul (2) 
^ k 

donde son los autovalores del espectro clasico de energi'as, correspondientes a autofun- 
ciones de cuadrado integrable, y la degeneracion de cada uno de ellos. 

La suma en (|^) resulta, en general, divergente. A fin de dar una interpretacion a las 
divergencias y poder calcular (0) debe recurrirse a alguna tecnica de regularizacion. En este 
trabajo, usaremos el metodo de regularizacion de la funcion zeta |^,|[. 

La funcion zeta de un operador A es definida, cuando son conocidos sus autovalores A„, 
como 

a = EV- (3) 

n 

Esta suma resulta convergente solo para Re{s) mayor que cierta cantidad, que depende 
de la dimension del espacio y del orden del operador 0. Sin embargo, la funcion puede 
extenderse analfticamente a todo el piano complejo, como una funcion meromorfa con polos 
en un conjunto discreto de valores de s. 

En nuestro caso, definiremos la suma @ como 

^0= (-^E^^i^^r^) , (4) 

donde /i es un parametro de escala con unidades de energia. 

Es decir: con el objeto de evaluar las energfas de vacfo que aparecen en la definicion (|1|) de 
la energia de Casimir, realizaremos la extension analftica de la funcion zeta correspondiente 
al hamiltoniano clasico, evaluandola posteriormente en el valor s = — 1 de su extension 
analftica. Al hacerlo, pueden ocurrir dos cosas: que el procedimiento arroje un resultado 
finito o bien que la extension presente un polo en el valor de s de interes. Cuando esto ultimo 
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ocurre, debe realizarse una renormalizacion de los parametros en el hamiltoniano clasico, a 
fin de dar una interpretacion fisica al resultado. En lo que sigue, veremos ejemplos de ambas 
situaciones. 

El contenido del presente trabajo es el siguiente: 

En la seccion II, determinamos el operador hamiltoniano clasico para fermiones en un 
campo magnetico uniforme en la direccion del eje de acuerdo al gauge elegido para el 
potencial electromagnetico y la representacion elegida de la matrices de Dirac. 

En la seccion III, obtenemos los autovalores del espectro clasico de energia, en caso en el 
que los fermiones estan confinados en un piano perpendicular al campo magnetico, hallando 
las autofunciones de cuadrado integrable del hamiltoniano clasico con impulso = 0. 

En la seccion IV, obtenemos las energias de vacio y de Casimir del problema en 2 + 1 
dimensiones, mediante la tecnica de regularizacion de la funcion zeta. 

En la seccion V, determinamos el espectro clasico de energia cuando todo el espacio 
tridimensional es accesible para los fermiones. Lo hacemos aplicando una transformacion de 
Lorentz sobre las autofunciones obtenidas en la seccion III. 

En la seccion VI, obtenemos las energias de vacio y de Casimir del problema en 3 + 1 di- 
mensiones, por medio de una regularizacion zeta y una renormalizacion de la carga electrica. 

En la seccion VII, consideraremos la variacion en la energia de vacio por efecto conjunto 
del campo magnetico y condiciones de contorno en la direccion del campo. Aqui nuevamente 
relizamos una regularizacion zeta y una renormalizacion de la carga electrica. 
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II. GENERALIDADES 



Estudiaremos el problema de un campo de Dirac en interaccion con un campo electro- 
magnetico, caracterizado por su potencial Af^{x). La ecuacion de Dirac en el espacio de 
Minkowski de 3 + 1 dimensiones y en unidades naturales, h = c = 1, se escribe entonces 



ih'^d^ - ^l^^ix - m)* = 0. 



(5) 



Trataremos el caso de un campo magnetico uniforme B en la direccion del eje z. Sabemos 
que la eleccion de un potencial no es unica y que el operador hamiltoniano depende de esta 
eleccion. Eligiendo el potencial electromagnetico = (0, — fx, 0) los operadores hamil- 
toniano y componente z del momento angular total conmutan y, por lo tanto, tendremos 
autofunciones simultaneas de ambos operadores. 



Utilizaremos la siguiente representacion de las matrices de Dirac. 



U o] 


( 




[ a,) 


V 



-i(T2 
— i(72 

donde cr, son las matrices de Pauli. 



^ i(Ti ^ 



U tax 



7 = 



iai 



(6) 



A fin de determinar el espectro clasico de energias, factorizamos la dependencia temporal 
del espinor de Dirac y lo escribimos en terminos de espinores de dos componentes 



-iEt 



( 0(f) ^ 



-iEt 



(7) 



Obtenemos entonces 



H 



(j){x] 



(8) 



donde el operador hamiltoniano esta dado, en el caso de campo electrico nulo, por 



H = 7°7(-iV - eA) + 



(9) 



con 



m 

eB, 



id^±dy- ^{y±ix 
o, en coordenadas cilmdricas 

{dr ± ± f r) 



H. 



± 



±ie 



—m 



(10) 



(11) 



En lo siguiente estudiaremos, en primer lugar, el problema de un campo de Dirac confi- 
nado en un piano perpendicular al campo magnetico. A continuacion analizaremos el caso 
en el que todo el espacio es accesible a los fermiones y, finalmente, estudiaremos el efecto de 
condiciones de contorno en la direccion del campo. 
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III. ESPECTRO EN 2 + 1 DIMENSIONES 



Si los electrones estan confinados en un piano perpendicular al campo magnetico, en- 
tonces las autofunciones resultan independientes de z 

\ 



lEt 



Xir,9) 

y el operador hamiltoniano resulta diagonal en bloques 



V 



(12) 



H 



H+ 
H 



J 



(13) 



con H± igual que en la ecuacion (O). 



Cada uno de estos bloques corresponde a una de las representaciones no equivalentes de 
las matrices de Dirac en 2 + 1 dimensiones [Bf . 



Estudiaremos, en primer lugar, el problema 



m 



ie-'^ (dr - Ide - ^ 



le 



id 



-m 



{r,e) = E(t){r,e). (14) 



Dado que vamos a obtener una base comiin de autofunciones de los operadores hamil- 
toniano y componente z del momento angular total, proponemos a autofuncion de este 
ultimo 



V 



, con A; G Z. 



Entonces, se satisface el siguiente sistema de ecuaciones de primer orden acoplado 

^ d k + 1 eB \ ^ r / \ 

■ -rrr gk[r) = [E - m)tk[r) 



(15) 



(16) 



z I - - + I /fc(r) = {E + m)gk{r), 



^dr r ' 2 ^ 
o bien, definiendo la variable adimensional p = 

7 



(17) 



1/2 



p]9k{p) = {E-m)fk{p) 



(18) 



(19) 



Analizaremos, en primer lugar, el problema E ^ —m. En este caso deberemos resolver 



pdp\dp) p2 P 



fkip) = 



(- 



E^-m^ + eB{k + l) Mp), 



(20) 



9k{p) 



(21) 



E + m \ 2 / ydp p 
A fin de resolver la ecuacion de segundo orden, hacemos el siguiente cambio de variables 



y proponemos soluciones de la forma 



Id ^d d ^^d 

pdp ^ ^ dp d^' 



(22) 



fk{p{0) = e-¥F{0. 



(23) 



Al reemplazarlas en la ecuacion diferencial, obtenemos 
2 



eB 



E^ - + eB{k + 1) e"5^2F(^) 



rl f d\ A;2 ; 



/ d'^ d k'^ \ ^ k , , , , 



Dado que 



I (e-dF(O) = e- 



" 2 ^ 2 ^ dr^^ ' ' 



(25) 



de 



rf^F fk ' 



dF (l 1(1- 
l4 2^^ e 



(26) 
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F{^) satisface la ecuacion diferencial 



2 



que puede identificarse con la ecuacion de Kummer 



donde b= k + ly a = 

Las soluciones de esta ecuacion, estudiadas detalladamente en el Apendice A, conducen a 
autofunciones de cuadrado integrable solo en los casos en que a es un entero negativo o cero, 
J b > a + 1. Estas soluciones son polinomios generalizados de Laguerre. Cuando a = —n 
los autovalores permitidos para la energfa deben satisfacer 

E^ = m'^ + 2eBn con n = 0,l,2,... (29) 

mientras que las soluciones F{C,) resultan 

F{0 = { " , (30) 

[ Ck,n.r'L;^1kiO si~n<k<0 

donde Ck,n es una constante compleja que se ajusta en cada caso para que los estados resulten 
de norma 1. Entonces, las soluciones de la ecuacion (EOl) son 



fk{p) 



Ck,n-e--p-''L-'li^{p'^) si-n<k<0 



Usando relaciones de recurrencia entre los polinomios generalizados de Laguerre es 
facil obtener las soluciones de la ecuacion (|2ll) , para n 

9k{p) = S , 1 2 > (32) 

[ ^C,,J^n.e-^p-^-^L-l-,\p^) si -n < k < 

y 9kip) = para n = 0. 
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Ahora analizaremos el caso E = —m. En este caso, la ecuacion ([TOD queda desacoplada 

( d _k 
[dp p 



^-- + p]fkip) = 0, (33) 



y su solucion es simplemente 

fkip) = C,.e-^p\ (34) 

donde k solamente puede tomar valores enteros positives para que la componente superior 
de (j) resulte de cuadrado integrable, y Ck es una constante de normalizacion. 
Por su parte, la ecuacion (|I8|) es una ecuacion de primer orden inhomogenea 

f d k + 1 \ . . 2miCk _Z k /or^ 

Proponemos soluciones de la forma 

9k{p) = G(p)e^p-('=+i), (36) 

por lo que G{p) debe satisfacer 

dG _ 2mtCk 2k+i /o7^ 
dp (—) 



G{P) = / dpe-^'p''^\ (38) 

V 2 ) 



y haciendo el cambio de variables p^ = x 



G{p{x)) = ^ / dxe "a; 



miC, ^ + Cte. (39) 



Por lo tanto, la solucion general de (p5|) es la solucion general de su homogenea corres- 
pondiente mas una solucion particular de la inhomogenea 



miCk A k\ 



10 



la cual conduce, para todo valor de de la constante D, a autofunciones de la componente 
inferior de (p que no son de cuadrado integrable. 

Por lo tanto, los autovalores de que corresponden a autofunciones de cuadrado inte- 
grable son los siguientes: 

{m 
^ (41) 
±Vm2 + 2eBn con n = 1,2,3,... 

El problema 



H-xir, 



-r 

2 



m ie*^ {dr + ^dg + 

y yur - Ide - ^r) -m j 



x{r,9) = Ex{r,9), (42) 



se resuelve en forma enteramente analoga. 

En efecto, factorizando la dependencia angular, obtenemos un sistema de ecuaciones de 
primer orden acopladas similar al de (|1^) y (|1^), que puede desacoplarse en una ecuacion 



de segundo orden para el caso E ^ m. Nuevamente, los autovalores de la energfa satisfacen 

E^ = m^ + 2eBn con n = 0,1,2,... (43) 

y las autofunciones pueden expresarse en terminos de polinomios generalizados de Laguerre. 

Por ultimo, al estudiar las soluciones del caso E = m, se encuentra que este autovalor 
no tiene asociadas autofunciones de cuadrado integrable. 

Finalmente, nos quedan los siguientes autovalores de correspondientes a autofun- 
ciones de cuadrado integrable: 

-m 

E_ = { (44) 

±Vm^ + 2eBn con n = l,2,3,... 
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IV. ENERGIA DE VACIO EN 2 + 1 DIMENSIONES 



Como se explico en la seccion I, la energia de vaci'o esta dada por la ecuacion 

1 



^9n-\En\ 



(45) 



donde los En son los autovalores dados por las ecuaciones (^) y (|^). El factor de dege- 
neracion Qn proviene de los distintos valores que puede tomar k (cualquier entero mayor 
o igual que n). Este puede calcularse contando el numero de estados que se encuentran 
en el problema de la particula libre con modulo cuadrado del impulso lineal entre 2eBn 
y 2eB{n + 1), si pensamos que al encender el campo todos estos niveles colapsan en uno. 
Integrando la densidad de estados de la particula libre, resulta 



9n 



LxLy 



J J dp^dpy 



(27r 



A ^ r^2eB{n+l) A „ 

iTT / pap = —eB, 

zvr 



V2eBn 



(46) 



donde A es el area perpendicular al campo magnetico. 

Conviene, entonces, considerar la energi'a de vacio por unidad de area 

^0 e5^, , eB f ^ ^ \ 

£o = = --r E \En\ = -—2 m + 2 ^ Vm2 + 2eBn . 
A 47r „ 47r V J 



(47) 



Como es usual, este resultado es solo formal y para darle sentido debemos regularizarlo. 
Hacemos esto con el metodo de regularizacion de la funcion zeta, descrito en la seccion I. 
Debemos multiplicar a la suma por un parametro para que conserve unidades de energia 
por unidad de area. Usamos /i con unidades de masa 



2eB 

4:71 



i + 1 



m 



+ 2 ^ (m^ + 2eBn 



n=l 



2eB 
~2^ 



m 



n=0 



-s/2 



m 



2eB 



s/2 



(48) 



s=-l 



Es conveniente introducir los parametros adimensionales 



b = 


2eB 






c = 


9 
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(49) 



de mo do que 



. ^ n=0 ^ " 



(50) 



s=-l 



Para Re{s) > 2, la serie en la ecuacion (^) converge y coincide con la funcion zeta 
Hurwitz. 



Cnis^a) = J2{n + ay 



(51) 



n=0 



que es convergente para Re{s) > 1, pero se extiende analiticamente a una funcion meromorfa 
con un polo simple en s = 1. Por lo tanto, la densidad de energia de vacio puede definirse, 
por extension analitica, como 



--(rn«(.v2.^ 



(52) 



S=-l 



Dado que (h (■5/2, |^ tiene un linico polo simple en s = 2, la densidad de energia de 
vacio por unidad de area resulta finita e igual a 

2eB 



27T 



- - {2eBf'\H -1/2 



2eB 



(53) 



Vale la pena sefialar que, en el caso de un campo puramente magnetico, la energia de vacio 
debe coincidir con la densidad lagrangiana efectiva al orden 1 loop del sistema. En efecto, 
el resultado de la ecuacion (|53|) coincide con la densidad lagrangiana efectiva obtenida en la 



referenda [jTO[, a partir del calculo con regularizacion zeta del determinante del operador de 
Dirac. 

Para calcular la densidad de energia de Casimir, debemos conocer la densidad de energia 
de vacio del problema sin campo. A fin de evaluar (^) para 5 = 0, usaremos un desarroUo 
asintotico de la funcion C//(s,a), valido cuando |a| — » oo |TT . 

1 



CH{.s,a) 



Vis) 



a-(«-i)r(s-l) + -a-T(s) + 



+ E 



2n 



-(2n+s-i)p^2n + s - 1) + 0(a-('"+^+^)) 



(54) 
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donde Bor, son los coeficientes de Bernoulli. 



Al escribir la densidad de energia de vaci'o en terminos de este desarrollo asintotico 
result a, para eB <^ rn? 



m m 



1 



^0 



E 



Bo„ /2eB\ 



2n 



j T{2n - 3/2) + O 



/2e5\2("+^)' 



(55) 



Stt 2tt r(-l/2) {2n)\ \ irfi 

3 

Vemos que el primer termino, independiente del campo magnetico, es y por lo tanto 
corresponde a la densidad de energia de vacfo del problema sin campo. En el Apendice B 
obtenemos el mismo resultado sumando directamente los autovalores del espectro clasico de 
energi'as de la particula libre. 



Resumiendo, la densidad de energia de Casimir resulta, para valores arbitrarios del campo 



^Cas — ^0 — £^0 



2eB 
~2^ 



--(2eB)'/^C« -1/2.^ 



3^' 



(56) 



Notese que, para campos pequefios, el segundo termino del desarrollo ( |53D cancela, en la 
ecuacion (|53|), el termino proporcional al campo. De este modo, resulta para la energia de 
Casimir un desarrollo en potencias pares del campo magnetico. 
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V. ESPECTRO EN 3 + 1 DIMENSIONES 



En lugar de resolver explicitamente (Q) podemos aplicar una transformacion de Lorentz 
sobre las soluciones independientes de z. En particular, considaremos un boost en la di- 
reccion del eje desde el sistema de referenda en el que = hasta el sistema en el cual 



^ cosh C(j sinh uj ^ 



10 
1 
sinhtt; coshcj 



donde cosh uj 



Vi 



Entonces 



E' = coshtuii^ 
p'^ = sinhiX'E' 



por lo tanto 



E'^ = cosh' uE'' = (1 + sinh' uj)E'' = E'' + p'/. 



,2 . , 772 



2,^/2 



Por su parte, los espinores se transforman segiin 



donde 



S = exp (-.| ^[7°, 7^]) = cosh (|) U + sinh (|) 7%^ 
En la represent acion elegida para las matrices de Dirac resulta 

S 



^ cosh(f)l2 -sinh (1)^2^ 



A partir de las soluciones 



y - sinh (f ) cxa cosh (f ) I2 



<f){r,e) 




-lE+t 



(57) 



(58) 



(59) 



(60) 



(61) 



(62) 



(63) 
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donde es una de las autofunciones del bloque se obtiene 



/ 



cosh (f ) 0(r, 9) 



\ 



V 



-iE+t 



sinh (f ) cT20(r, 6') 
o bien, en terminos de las coordenados del sistema primado, 



(64) 



cosh (f ) 0(r', ^ 
■sinh(f)(720(r',^')y 
Estas autofunciones corresponden a los autovalores de la energia 



(65) 



E' 



con p'j. e R y n = l,2,3, ... 



(66) 



Similarmente, a partir de las autofunciones x del bloque se obtienen los siguientes 
espinores en el sistema de referenda primado 



^-sinh (f) a,x{r',e')^ 



\ 



cosh (f)x(r',^') 



(67) 



que corresponden a autovalores 



E' = <^ 



/2 



con p'^ e R y n = l,2,3, ... 



(68) 
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VI. ENERGIA DE VACIO EN 3 + 1 DIMENSIONES 



En este caso, la energia de vacio se obtiene a partir de sumar los autovalores dados en 



la ecuaciones (|66D y 



(69) 



donde el factor de degeneracion g es nuevamente -^eB. 



Redefinimos el mdice de la suma para que esta comience desde y definimos la densidad 
de energia de vacio por unidad de volumen 

Eq eB 



r dp, (\/m'^+pl - 2 ^m"^ + pI + 2eBn \ . (70) 

V n=0 / 



AL, (27r)2 J-c 

Nuevamente, la expresion que encontramos para la densidad de energia de vacio es formal- 
mente divergente. Repitiendo el procedimiento antes utilizado, la extendemos analiticamente 
con la tecnica de regularizacion zeta. Como antes debemos agregar algiin parametro para que 
la densidad de energia de vacio mantenga las unidades apropiadas. Usando con unidades 
de masa, resulta 

^ \^^'^' / _ dp, ^m^+pi)~"^-2 

n=Q J J s=-l 



/i^+i r dp, (m^ + p2)-«/2 _ 2 £ (m^ + p2 + 2eBn 



-s/2 



(71) 



Si usamos la transformada de Mellin 12 



T(s) Jo 



s-l-Zt 



dtf-'e 



(72) 



la ecuacion ( |7TD se escribe 



^0 



eB 
(2^ 



/oo 1 poo 

^P^FTT^/ dtt^/'-' exp[-{m' + pl)t] - 
-oo 1 {s/2) Jo 



1 /.oo , 1 "I 

-^EwTa dtt'^'-'eM-irn'+Pl + ^eBn)t] . (73) 

n=0 r(s/2) Jo J J 

El orden de integracion puede intercambiarse para Re{s) suficientemente grande, de 
mo do que 
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So 



eB 



s+1 



3/2) JO J~oo 



T{s/2) 



°° 1 I 
-'^T.T^TTa dtf^'-'exp[-{m^ + 2eBn)t] dp,expi-plt) \ 

„=o r(s/2) 7o j-oo J J 



Usando el resultado de la siguiente integral definida 



dxe~ 



resolvemos las integrales a lo largo de pz- Asi obtenemos 



eB 



fx 



s+l 



r(s/2) Jo 



dtt 2 ^exp(-m t) 



(74) 



(75) 



-2E 



^,Tis/2) Jo 
Luego de realizar las integrales en t, result a 

eB 



dtt^~'^ exp \-{m^ + 2eBn)t 



£o 



(2vr) 



p/ g— 1 N OO 

2^ii-l-(^(m2 + 2e5n)- 



s-l 
2 



(76) 



(77) 



s=~l 



r(V2)^ ^ r(./2)^„^ 

Nuevamente, la serie en el segundo termino define, para Re{s) > 3, una zeta de Hurwitz. 
Por extension analitica, obtenemos 



2eB 



TT < /i 



rf— 1 



r(./2) 



, s-i /s — 1 \ 



(78) 



s=-l 



A diferencia de lo que ocurria en 2 + 1 dimensiones, la extension de la funcion zeta 
del hamiltoniano clasico presenta, en 3 + 1, un polo simple en el valor s = — 1, debido al 
comportamiento de la funcion r(s) en los enteros negativos. 

Usando nuevamente los parametros adimensionales definidos en (|^), ( [78D puede es- 
cribirse 



£o 



{2tt 



s-l 
C 2 



r(^; 



T{s/2){s + l){s-l) 



,_£^^ fS-1 C 



3-1 

C 2 s-l 
2 ( 



s=-l 



-1 C 

"2~'6 



(79) 



s=-l 



o bien 



TT^ S + 1 



s=-l 



donde 



C 2 ,_£^ . /s — 1 c 



r(./2)(.-i) 

A fin de identificar las partes singular y finita de la densidad de energia de 
un desarroUo alrededor de s = — 1. 



1-^ 



z{-i) 

s + 1 



+ Z'(-1) 



s=-l 



Para la parte singular debemos calcular 

r(i) 



Z{-1) = 



-2)r(-i/2) 

y usando la extension analitica de la funcion Zeta de Hurwitz 



Cif(-l,a) = (^a^ - a + i 



Z{-1) 



80r {b"^ 6 r 



La parte finita esta dada por 
b 



160r 



(^lnc-21nKif(-l,0+2C;/(-l, 



^(^(l)-^(-l/2) + l)g + [ 



160r 



-1. 
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Reemplazando (|85|) y (|86D en (|8^), la densidad de energi'a de vaci'o resulta 



+ 



6 



7r2(s + l) 87r2(s + l) 



+ 



+ 



167r2 



(^(l)_^(_l/2) + l-lnc) + 



+ 



+ ic^ - cb+ — Uni- 



(87) 



Vale la pena sefialar, en este punto, que la parte finita de la expresion anterior coincide, 
a menos de terminos proporcionales a 6^ (y, como veremos, renormalizables) con la densidad 
lagrangiana efectiva calculada en |JlO|| a partir del determinante del operador de Dirac. La 
diferencia entre ambos resultados esta de acuerdo con la discusion general presentada en 



12 



De la expresion (p7|), vemos que el residuo en el polo presenta un termino independiente 
del campo magnetico, y otro proporcional a su cuadrado. 

A fin de dar una interpretacion fisica a este resultado verificaremos, en primer lugar, que 
el termino independiente de h se cancela con uno identico que aparece en la densidad de 
energi'a de vacio a i? = 0. 

Para hallar la densidad de energfa de vacio sin campo £:[), po demos considerar un desa- 
rrollo de asintotico de la derivada de la funcion zeta en el limite de los campos pequeiios, y 
quedarnos con los terminos del resultado final independientes del campo (en el Apendice C 
la reobtendremos directamente a partir de la suma de los autovalores del espectro clasico de 
la particula libre). 

El siguiente desarrollo para la derivada de la funcion zeta de Hurwitz es valido en el 



Hmite a oo 10 



m—l 



-2n 



Crjf— 1, a) = - fa^ — a H — \hia H h , 

' ^ 2 V 67 4 12 ^^{2n){2n + l){2n + 2) 



+ 0(a-'™). (88) 



Haciendo uso del mismo resulta, para h <^ c 
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_ 1 

2 62 ~ 6 



m— 1 



-2?i 



_2 _B2n+^ 

(2n)(2n+l)(2n + 2) 



2m' 



y, por lo tanto, 



^7r2 s + 1 



s=~l 



+ T^[^(l)-^(-l/2) + 3/2-lnc] 



De este modo, la densidad de energia de Casimir resulta 



(89) 



(90) 



^Ca 



^0 



6 



ill' 



is + 1] 



+ 



167r2 



1,2 2 

(^(l)_^(_l/2) + l-lnc)--|- + 



6 



+ [c' -cb+-]\n{^] -2b\'H(-l, 



(91) 



Esta expresion presenta aiin una divergencia, cuya dependencia funcional con el campo 
es igual a la que aparece en la densidad clasica de energia 



(92) 



2 8e2 

A fin de eliminar la singularidad que persiste en la densidad de energia de Casimir, 
renormalizamos la carga e. Para ello, consideramos a los terminos proporcionales al cuadrado 
del campo como la correccion cuantica de la densidad de energia electromagnetica clasica. 

De esta manera, definimos la carga renornalizada a traves de 



d{ecid + So) 



62=0 



8e| 



(93) 



Evaluamos esta derivada a partir de la ecuacion (79) y obtenemos 



d{e,i, + e,) _ 4/r(^) 



(2vr) 



r(./2) 



s — 3\, s+i / s — 1 c 



S — 1\ ,_s+3 ^ 



s + 1 c\ c "2 b ^ 
~ir'b) ^ 4 



(94) 



s=-l 
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donde hemos usado la siguiente propiedad de la funcion zeta de Hurwitz 

dCHis,a) 



da 



-sCni-s + l,a) 



(95) 



Podemos evaluar las zetas de Hurwitz en el limite 6^ — con el desarrollo ([5^). Final- 
mente, result a 



8e2 8 



4 ^-4irv 2 
yU -D2 C 2 



r(./2) 



(96) 



s=-l 



Dado que esta expresion es singular en el limite s — > — 1, igual que como hicimos antes, 
separamos las partes finita y singular mediante un desarrollo alrededor de s = — 1, con lo 
cual obtenemos 



d{ecid + £o) 



fe2=0 



y de aqui se tiene 



4 4 

8e2 967r2 



s + 1 



+ ^(1) -^/'(-l/2) -Inc 



s=-l 



(97) 



1 + 



127r2 



_|_ + ^(l)_^(_l/2)-lnc 



127r2 



s + 1 



+ -^/'(-l/2) -Inc 



(98) 



Entonces, el parametro que aparece en el hamiltoniano debe ser renormalizado segiin: 



4 + 



127r2 



s + 1 



+ ^(1) - V'(-l/2) -Inc 



+ 0{e%). 



Finalmente, para la densidad de energia al orden 1 loop: 



£Cas — £dd + ^0 ~ ^0 



(99) 



8el 167r2 



(100) 



y en terminos de los parametros originales 
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£Ca 



, 1 ({2eRBy 



2 ^ 167r2 \ 6 



In 



+ m 



\2eRB) 2 



:ioi) 



Notese que, una vez renormalizada la carga, la densidad de energia de Casimir es inde- 
pendiente del parametro de escala /i. 



23 



VII. ENERGIA DE VACIO EN 3 + 1 DIMENSIONES CON CONDICIONES DE 

CONTORNO 



En esta seccion, analizaremos el efecto de tamafio finite en la direccion del campo sobre 
la energia de vacio. Impondremos sobre los campos de Dirac condiciones de contorno no 
locales, del tipo 



i^{r,9,L) = e^Xr,^,0). 



(102) 



Vemos que esta condicion se reduce a la condicion de period! cidad para a = 0, mientras 
que de a = TT se obtiene el problema con condiciones de contorno antiperiodicas. 

Nuevamente tenemos autovectores simultaneos de la energia y del impulso en la direccion 
del campo, solo que ahora el espectro de autovalores pz pasa a ser discreto. 



Si e*^^^ = e*" =^ Pz = ^^yi^ donde k es un entero • 



(103) 



Entonces, en la expresion para la energia de vacio, debemos ahora sumar sobre el espectro 
discreto de autovalores Pz 



eBA 



271 



E 

k=—oo 



9 / 27rfc + a 



1/2 



+ 2E 



n=l 



m 



27rk + a 
L 



+ 2eBn 



eBA g 



271 



k=—oo 



2E 

n=0 



o / 27[k + a\ 
+ + 2eBn 



1/2 



m 



'27ck + a' 
L 



1/2' 



1/2' 



• (104) 



Una vez regularizada mediante la tecnica de la funcion zeta, ( p.04| ) se escribe 



eBA 
2n 



s+l 



E 2E 

k=—oo n=0 



9 / 27Tk + a\ 
+ ( 1 + 2eBn 



-s/2 



n / 2Tik + a 

" + — ^ 



-s/2" 



(105) 



Usamos la transformada de Mellin (|72D para luego poder hacer la suma en n 



24 



^0 7^ 1 ^ 1 ^ 



27r 



k=—oo kn= 



r(s/2) 



X / dtf/^-' 
Jo 



exp 



+ + 2eBn t 



1 

3/2 io ^ 



r(./2) 



L 



+ I I t 



L 



s=-l 



27r 



^^^,V{s/2) Jo 



a 



+ {-] + 2eBn t 



X 



X 2^ exp I — —tk —tk 



k=—oo 



L2 L2 



- — — / dtf/'^-^ exp 
r(s/2) io ^ 



+ I ? ) It 



^ f Ana , Air^ , 
X 2^ exp I — r^t/c —tk 

k=—oo 



L2 L2 



(106) 



s=-l 



Podemos expresar este resultado en terminos de la funcion theta de Jacobi, que queda 



definida por [0] 



e{x,y) = e^'^'^e-'^''^. 



(107) 



k=—oo 



Por lo tanto 



^ f Ana , 2at Ant 

E^-p(— '*^j=»(-^.7? 



(108) 



Considerando la siguiente propiedad de inversion de la funcion theta de Jacobi []T3 

1 



^(x,.)^^expp: 4^,1 , 

\ y vv y 



(109) 



result a 
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J2 

fc=— oo 



-tk -tk^ 



W / L 



L2 L2 



t \2txJ ^ V L2 / I 27r' 47rt, 



E e^^(iak-—e]. (110) 



A;=— oo 



Reemplazando en la expresion para la energia de vacio 



eBA / ^^^J^g^ 



27r ' 27T 



o^{s/2) Jo 



dtt^~^ exp [- (m^ + 2e5n) t 



X E e^'^^exp 



A:=— oo 



At 



1 /-oo °o / 

-— / c/tt^-^exp (-m^t) ^ e^'^^exp - 



L'k 



;iii) 



s=-l 



Nuevamente, calculamos la densidad de energia de vacio por unidad de volumen. Escri- 

biendo por separado el termino correspondiente a k = podemos identificar a la densidad 

de energia de vacio sin condiciones de contorno. 
T. ei? / „ , 1 ^/7^ 



(2^) 



- 2 E / dtt^~^ exp [- {n? + 2eBn) t 



n=0 



+ 



s=-l 



(2vr)2 V" r(./2) 



{/•oo oo / 

y dtt— -1 exp(-m2t)2 E e'''" exp ^ 



^ r°° -1 r / \ n °° / 

- 2 E / exp [- (m^ + 2eBnj tj 2 ^ e'''" exp - 

n=0''^ k=l \ 



At 



;ii2) 



Efectivamente, los dos primeros terminos corresponden a la densidad de energia de vacio 
£0 del problema en todo el espacio estudiado en la seccion previa (ver ecuacion ([76|)) 

2eB \ ^fi'+^ 



^0 ~ ^0 



(27r)2\ r(s/2) 1^,^ Jo 



Aka 



dtt 2 ^exp(-m t) X 



X 



oo 
. n=0 



2eBt\n 



At 



exp 



(113) 
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A fin de estudiar la variacion de energi'a debida al tamano finito, realizamos la suma 
geometrica, cuyo resultado es 



1 



n=0 



— X 



para a; < 1, 



(114) 



y evaluamos para s = — 1 



^0 — 



eB 



POO 

e'''" / dtt-'^ exp{-mH) coth{eBt) exp 
k=i •'^ 



At 



(115) 



Cabe senalar aqui que, para a = 0, vr, este resultado coincide con el encontrado, mediante 



el esquema de regularizacion del tiempo propio en |1^ . 

A fin de obtener una expresion mas explfcita para la densidad de energia de vacio, usamos 
el siguiente desarrollo para la cotangente hiperbolica 

„2 oo -j^ 



a;coth(a;) = 1 + 



X 



2x^y 



;ii6) 



La contribucion a ( |115| ) del primer termino del desarrollo es 



1 oo 
k=l 



(2vr)2 



^'^^ 1^ dtt-^'expi-mH-^ 



1 ifca /""^ . / 

- — > e / axxexp 

(27r)\4l Jo X 



L^k X 



;ii7) 



La integral puede resolverse con el siguiente resultado [0, el cual es valido para Re{(3) > 
y i?e(7) > 

u/2 



dxx''-^ exp 1^-^ - = 2 1^^^ K, (2 



(118) 



donde K^, son las funciones de Bessel modificadas. La contribucion de este termino, inde- 
pendiente del campo, es 



2^2 00 ^ika 



Y^—K,{mLk). 



;il9) 



Este termino se reduce, para a = 0, al resultado bien conocido para la energia de Casimir 



de un campo de Dirac en una caja con condiciones periodicas |jT5| 
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El segundo termino del desarrollo ( p.l6| ) dara una contribucion cuadratica en los campos 



3 (27r) 



Y.^^'^ \ t^tt"^ exp j^-m^t 



4t 



3t)|:^-A«(-««). 



(120) 



donde nuevamente hemos usado (|118|) para resolver la integral. 

Podemos analizar la convergencia de estos dos primeros terminos a partir del compor- 



tamiento asintotico de las funciones de Bessel modificadas para argumentos grandes [jT4 
Para z ^ oo 



l + 0{- 

z 



;i2i) 



independientemente de u. Vemos que el factor e ^ garantiza la convergencia de las sumas 

Oy dUD 



Para evaluar el ultimo termino en (|116|) sumamos y restamos el termino que correspon- 
deria a = para poder extender la suma desde — oo hasta +oo y poder reconstruir la 
funcion theta de Euler. 



,{eB)_ 
' (2vr)2 



4 oo 



fc=i 



ika 



exp (^—m?t 
'^^''^S7r2g2[(e5t)2 + 7r2g2 



4t 



\4 oo 



(27r)2 



fc=— oo 



(■oo °° 



exp ( — m^t 



4t 



>V[(e5t)2 + 7r2g2 



+ 



+ 



{eBY 



exp^—mH) 



122) 



(27r)2 Jo TT^q^[{eBt)^ + n^q^] ' 

En el primer termino podemos formar el siguiente factor y reescribirlo usando la definicion 
de la funcion 9{x,y) y su propiedad de inversion ( |109| ) 



Aka , 



fc=— oo 



t (In 



71 \L 



E exp 



A:=— oo 



2k7i + a 
L 



(123) 
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La contribucion a la densidad de energfa de vacio es 



{eB) 



2 oo oo 

E E 



1 



fc=— oo <J=1 ^ 



exp 


— m^t — 




2 " 






2' 





+ 



+ 



2 oo 

E 



1 f°° ^f-^ exp(— m^t) 



27r / ^ 7r2g2 Jq 



(124) 



Las integrales las resolvemos usando la siguiente solucion que es valida para Re{(3) > 0, 



ReiiJ,) > y i?e(z/) > -1 [T] 



x'^exp(— /ix) r(z/) 



X2 + /52 



exp ( + i- — ^ I r(l — u, ifij3) + 



+ exp \—iix(3 — i- — ^ ] r(l — u, —ifijS) 



(125) 



donde T{a, z) es la funcion gamma incompleta 



La contribucion a la densidad de energfa es 



;5)3/2 g g 1 



87r5/2L . ^ ^ g3/2 



exp 



iq-K 



r, / 2A;7r + a 



27r 
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X F -1/2, 



r, / 2A;7r + a 



+ exp 



2A;7r + a 
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T 



r -1/2, 



iqT{ 
"eB 



9 / 2A;7r + a 



+ 



, {eBf ^ 1 



<?=! 



/ iiv?qTi\ ( iiv?qTi\ 
exp — r 0, — + exp 



eB 



eB 



^ *r 0, ^ 



eB 



eB 



(126) 



Para verificar la convergencia de las series, consideramos el comportamiento asintotico 



de las funciones gamma incompletas Para 2 — > oo en | arg(2;)| < 



1 + 



a-1 (a-l)(a-2) 



(127) 
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Para kj/oq suficientemente grandes, los sumandos se comportan de la siguiente manera: 
en ambos casos los primeros terminos del comportamiento asintotico se cancelan entre los 
dos terminos de cada sumando. Para el sumando de la primera suma (doble en este caso) 



oi7r/4 



(eB) 



5/2 



i4 ■ 



_|_ ^2fc7r+ay 

Para el segundo sumando, que solo se suma en q 

{eBf 



5/2 ■ 



(128) 



^2^2g4 



(129) 



Concluimos que las sumas de estos dos liltimos terminos de la densidad de energi'a de 
vacio son convergentes. La expresion que finalmente damos para la densidad de energi'a de 
vacio del campo de fermiones en interaccion con un campo magnetico uniforme y con una 
condicion de contorno en la direccion del campo es la siguiente 

= eo + e^, + - j E e^'°^o(mLfc) - 



87r5/2L ,,f-_, g3/2 1 \ e5 



k= — oo q=l 



X r -1/2, 
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iqn 
'eB 



9 / 2k7i + a 
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iqn 
"eB 
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2k'K + a 
L 



+ 



( iw?qT{\ ( irn?qTi\ 
exp — r 0, — + exp 



im?qTi' 



r 0, 



im?qTT^ 



. (130) 



eB J \ eB J "V eB ) \ ' eB 

En esta expresion de la densidad de energfa de vacio hemos aislado una contribucion, Eq, 
independiente de las condiciones de contorno y otra, Sq, debida exclusivamente al tamaiio 
finito e independiente del campo. Si definimos la densidad de energfa de Casimir como la 
diferencia entre Sq en ( p.30|) y su valor para 5 = y L — > oo,obtenemos (usando Eq de ([87|)) 
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eB 



Nuevamente, los terminos proporcionales a B^ pueden eliminarse por renormalizacion de 
la carga electrica. En este caso, debido al tamafio finito, apareceran nuevas contribuciones 
a la carga renormalizada (por ejemplo, el tercer termino en ( p. 301 )). Para ver que, en efecto, 
este es el linico termino que agrega una nueva contribucion a la carga renormalizada, la 
calculamos expHcitamente a partir de su definicion (pi) y de (|115|). Para ello, debemos 



mostrar que la serie derivada en ( 115 ) es uniformemente convergente. El segundo termino 
en ( |115| ) se escribe, al desarroUar la cotangente hiperbolica. 
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(132) 
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La serie derivada result a 

2e2 °° 



poo 



exp — m t — 



(133) 



(27r)2 Jo [(e5t)2 + 7r2g2]2 ^ ■ - " 4^ ^ 

y el modulo de sus terminos puede acotarse con los terminos independientes de B de una 
serie convergente 
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Concluimos que la serie derivada en ( |115|) es uniformemente convergente y, por lo tanto, 
podemos evaluar su derivada termino a termino. Asf, obtenemos 
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k=l 



De esta manera, se agrega un termino en la carga renormalizada que depende del tamano 
finito en la direccion del campo 
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La expresion que finalmente nos queda para la densidad de energia de Casimir 
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Otra vez debe destacarse que, una vez renormalizada la carga, la densidad de energia de 
Casimir no depende del parametro de escala /i. 



■ (137) 
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VIII. RESUMEN DE LOS RESULTADOS OBTENIDOS 



En este trabajo, hemos calculado las correcciones cuanticas a la energi'a de vacio (ener- 
gias de Casimir) para campos de Dirac en presencia de un campo magnetico uniforme de 
background B. Hemos analizado las correcciones debidas exclusivamente a la presencia del 
campo de background y al efecto combinado de dicha presencia y la imposicion de condiciones 
de contorno. En todos los casos hemos utilizado el metodo de suma de modos, combinado 
con la regularizacion zeta 

En primer lugar, hemos evaluado la energfa de Casimir para campos de Dirac confinados 
en un piano perpendicular a la direccion de B (problema en 2 + 1 dimensiones) . En este 
caso, la energia de Casimir, definida por extension analftica, resulta finita y no requiere, por 
lo tanto, una renormalizacion de los parametros de la teoria. 

A continuacion, hemos tratado el problema en 3 + 1 dimensiones, cuando todo el espacio 
tridimensional es accesible a los fermiones. En este caso, la correccion cuantica a la energi'a 
de vacio resulta divergente, con residuos en los polos de dos tipos: independientes de i? y 
proporcionales a B^. La primera de estas divergencias desaparece en el calculo de la energi'a 
de Casimir, dado que esta presente tambien en ausencia del campo magnetico. Los terminos 
proporcionales a B^ pueden eliminarse por renormalizacion de la carga electrica. 

Finalmente, hemos tratado el efecto combinado de campo de background y condiciones de 
contorno no locales ("twisted") [Q, cuando el espacio accesible a los fermiones tiene tamaiio 
finito en la direccion de B. En esta situacion, hemos podido separar las contribuciones a la 
energi'a de vaci'o debidas a uno y otro factor externo exclusivamente de aquella debida a la 
accion combinada de ambos. Aquf, nuevamente, debimos renormalizar la carga electrica a 
fin de obtener un resultado finito para la energfa de Casimir 

En los casos sin borde, hemos verificado la consistencia con resultados previos de otros 



autores [T^ para la accion efectiva (ya que para campo electrico nulo debe coincidir con la 



energi'a de vacio). En el caso con condiciones "twisted" hemos mostrado que, en los limites 
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particulares de condiciones periodicas y antiperiodicas, nuestra energfa de vacio, expresada 
en forma integral, coincide con resultados anteriores ||I6|. En este trabajo, hemos encontrado, 
para un parametro de "periodicidad" arbitrario, una expresion mas explicita de la energfa 
de Casimir, como serie de funciones especiales. 
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APENDICE A: SOLUCIONES DE LA ECUACION DE KUMMER 



Las soluciones independientes de la ecuacion de Kummer 



(Pw ,dw ^ 



se pueden expresar en terminos de 

, , , az a(a + l)z'^ a(a + l)...(a + n — l)z'^ 

definida solo para b distinto de un entero negativo o cero y 



(138) 



(139) 



U{a, b, z) 
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sinirb 



M{a,b, z) 



— z 



i_(,M(l + a-6,2-6,z) 



(140) 



_r(l + a-6)r(6) r(a)r(2-6) 

si bien estas dos funciones son soluciones de la ecuacion de Kummer, no siempre resultan 
linealmente independientes. 

Para estudiar los distintos casos consideraremos las cuatro siguientes soluciones 

yi{a,b,z) = M{a,b,z) 
y2{a, b, z) = z^-^M{l + a-b,2-b,z) 
ysia, b,z) = U{a,b,z) 
y4{a, b, z) — e^U {b — a, b, —z) 

Estudiaremos los comportamientos asintoticos de las soluciones en los distintos casos, 
para 2; — > oo y para 2; — > 0, y veremos para que valores de a y 6 conducen a autofunciones 
de cuadrado integrable, recordando que en nuestro problema b solo toma valores enteros. 
En particular consideraremos el comportamiento de M(a, 6, z) para z grande, cuando a es 
distinto de un entero negativo. 



(141) 



M(a, 6, z) = ^^e^z^-' 
r(a) 



1 + 0(|2;| ) para z ^ 00, Re{z) > Q 



(142) 



y el comportamiento de C/(a, 6, z) en un entorno del origen 
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U (a, b, z) 



l--z-^ + 0(\\nz\) 



T{a) 



si 6 > 2 



si 6 = 2 



Yj^[\\iz + %l){a)] + 0{\z\Y).z\) si 6=1 



^ + 0{\z\riz\) si 6 = 



r(l+a) 



r(i-b) 

T(l+a-b) 



+ 0{\z\) si6<0 



(143) 



1. a 7^ —n (negativo o cero) y b entero b> 1 
La solucion general es 



w{z) = Ayi{a, b, z) + 57/3(0, 6, z), 



(144) 



ya que su wronskiano es no nulo 



r a) 



(145) 



Si z — > 00, entonces yi ~ e^z"' ^ y, por lo tanto, las autofunciones que se obtienen a 
partir de esta solucion divergen exponencialmente. 

Si 2: — *■ 0, 1/3 ~ z^^^ para 6 > 2 y entonces las autofunciones, que divergen como 
p^^ para A; > 1 en el origen, no son de cuadrado integrable. Mientras que para 
6 = 1 la divergencia es logarftmica y en consecuencia esta componente del espinor si 
es de cuadrado integrable. Sin embargo, al derivar la componente restante de (|TU]) la 
divergencia en el origen va como k~^. 



2. a 7^ — n (negativo cero) y b entero b < 
En este caso, la solucion general es 
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w{z) = Ay2{a, b, z) + Bys{a, b, z), (146) 

dado que son soluciones linealmente independientes 

'^fe.!'3) = -f^f^.-V (147) 

Si z — > oo, entonces 1/2 ~ e^^;""** y ,por lo tanto, las autofunciones que se obtienen a 
partir de esta solucion divergen exponencialmente. 

Si 2; — > 0, 7/2 1 y entonces las autofunciones divergen como p'' con A; < —1 en el 
origen y no son de cuadrado integrable. 

3. a — —n (negativo cero) y b entero b> 1 
Escribimos la solucion general como 

w{z) = Ayi{a, b, z) + By4{a, b, z), (148) 

para las cuales 

W{y,, y,) = ^^^f""'^''^' (149) 

Entonces yi = M[—n,b,z) = ^fj_^-^-^^^_^^' (fc_x+„) -^n~''^; polinomio de Laguerre general- 
izado que conduce a autofunciones de cuadrado integrable. 

Mientras que 1/4 = e^U{b + n,b, —z) ~ e^{—z)~''~"' diverge exponencialmente, para 
z 00. 

A. a — —n (negativo cero) y b entero —n + 1 < 6 < 
La solucion general es ahora 

w{z) = Ay2{a, b, z) + By^i^a, b, z), (150) 
y el wronskiano de estas soluciones es 
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W^(l/2,l/4)=-^|^^-V (151) 

Tenemos = z^-^M{-n + 1 - 6, 2 - 6, = z^-^ ^^_,^^^Xi):S-^+n-i) con las 

que resultan autofunciones de cuadrado integrable. 

Por el contrario, las soluciones ^4 = e^U{b + n, b, —z) ~ e^(— 2;)"''"" divergen exponen- 
cialmente,para 2; — > 00. 

5. a — —11 (negativo cero) y b entero b < —n 

Para la solucion general nuevamente proponemos la combinacion 

w{z) = Ay2{a, b, z) + By^{a, b, z), (152) 
ya que wronskiano de estas soluciones es 

Si 2; — > cxo, entonces y2 ~ e^z""'^ y, por lo tanto, las autofunciones que se obtienen a 
partir de esta solucion divergen exponencialmente. 

Si 2; — >■ 0, 1/2 ~ 1 y entonces las autofunciones divergen como con k < —1 en el 
origen y no son de cuadrado integrable. 
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APENDICE B: CAMPO LIBRE EN 2 + 1 DIMENSIONES 



La energia de vaci'o del campo libre se obtiene de integral las energias de la parti'cula 
libre en el espacio de impulsos, multiplicadas por la densidad de estados de la particula libre. 
Esto es 



J 



Eo 1 1 

°" A " 2(27r)2 



P 



1 

2\ 2 



(154) 



Emplearemos el metodo de la funcion zeta para regularizar esta integral, el mismo que 
utilizamos en la seccion |IV|. El parametro /i con unidades de masa asegura que la expresion 
tenga unidades de densidad de energia por unidad de area. 

2 



^ (27r)2 
Usamos la transformada de Mellin (|72D 



/oo roo 
dpx / dpy{m'^ + pI + pI 
-oo J —oo 



2\-s/2 



s=-l 



(2vr) 



/i 



s+l 



T{s/2) 



Jo J—oo J~oo 



A s=-l 



y evaluando las integrales de los impulsos con 



(27r) 



r(s/2) Jo 



m t 



(155) 



(156) 



1 



(27r) 



s=-l 



(157) 



r(./2) 

Esta expresion para la densidad de energia de vacio esta regularizada ya que la funcion 
gamma solo presenta singularidades en los enteros negativos y cero. 



1 r(-3/2) 



m 



-m 



° (27r)r(-l/2)''" 37r 
y el resultado coincide con el termino independiente del campo del desarrollo (|55|) . 



(158) 
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APENDICE C: CAMPO LIBRE EN 3 + 1 DIMENSIONES 



De igual modo, expresamos la energia de vacio del problema de campo libre en 3 + 1 
dimensiones, 



En 



LxLyLz 



-2 • 2- 



1 1 



2 I „2\2 



d p Im + p 



(159) 



2 (27r)3 

Esta expresion es divergente. La regularizamos nuevamente por medio del metodo de la 
funcion zeta, e introducimos un parametro para que conserve las unidades apropiadas 

2 



(2vr)= 



/CO roo roo 

dpx / dpy / dp,{m^ +pI+pI+ V^^~"'^ 



(160) 



s=-l 



Hacemos uso de la transformada de Mellin (173) y obtenemos 



(27r)^ 



/i 



r(s/2) Jo 



X 



dp^e-P-^ / dpyC'Py^ / c/p^ 

OO J—oo J — oo 



-pit 



s=-l 



(161) 



Resolvemos las integrales con ( [75| ) y resulta 



(2vr) 



1 



/i 



s+l 



r(./2) 



7r^/2 / dtt^"^e 



s=-l 



4vr3/2 



(162) 



s=-l 



r(./2) 

Nuevamente usamos el parametro adimensional c definido por (^9]), y factorizamos la 
singularidad que la funcion r(^) tiene en s = —1, 



IT 



3/2 



r(^) 



r(s/2)(s-3)(s- l)(s + l) 



C 2 



(163) 



J s=-l 



y separamos en partes singular y finita haciendo un desarrollo alrededor de la singularidad 



TT 



3/2 



z{-i) 

s + l 



z'(-r 



(164) 



s=-l 



donde la funcion Z{s) queda definida por 
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Zis) 



r(^)c-^ 



r(./2)(s-3)(s-i)' 



(165) 



Para la parte singular 

Z{-1) 
Para el termino finite 

Z'i-l) = - 



r(-l/2)(-4)(-2) IQ^ 



(166) 
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V'(l)-^(-l/2) + --lnc 



(167) 



Finalmente, para la densidad de energia de vaci'o del campo libre en 3 + 1 dimensiones 
obtenemos 



m m 
^7r2(s + 1) ^ 16^ 



3 JTo2 

^(1) -^(-1/2) + --In — 

^ Li 



(168) 



que coincide con el resultado ( pOD obtenido en la seccion V. 
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